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1.0 SUMMARY 


The computational procedure and numerical results are presented for a hew method to 
solve Kiissner’s integral equation in the case of subsonic compressible flow about 
harmonically oscillating planar surfaces with controls. Since Kiissner’s equation is a 
linear transformation from- pressure to normalwash, the unknown pressure is expanded 
in terms of prescribed basis fimctions; the unknown basis function coefficients are 
determined in the usual manner by satisfying the given normalwash distribution, either 
coUocationally or in the complex least squares sense. The present method of solution 
differs fjom previoais ones im that the basis functions are defined in a continuous'' 
fasbion over a relatively small portion of the aerodynamic surface and are zero 
elsewhere. ,Tl)is method, termed the local basis function method, combines the 
smoothness and accuracy of di^ribution methods with the simplicity and versatility of 
panel methods. Predictions by the local basis function method for unsteady, flow are 
shown herein to be in excellent agreement with other methods. Also, potential 
improvements to the present method and extensions to more general classes of solutions 
are discussed. 


I 

k 

► 


i 


2.0 INTRODUCTION 


Aeroelastic analysis of flight structures involves the coupled dynamic response of 
structural and aerodynamical media. Within the confines of, lipear aeroelas^ci|;y, the 
coupling b determined by the matrix of energy that is transmitted pcross the bo4y-0u^ 
interface 'in the form of mechanical work done as the body deforms in one modp a^ainsf 
the pressure due to another mode. This matrix, called the generalized aerodynamic fqrcg 
’^coefficient' matrix, depends for its calculation upon knowledge of the unsteady pressure 
‘'diStributioh. 

• ‘i -t ^ - 


f 

i. 


i 


ft THte pft-eSsure 'Bi^ribution may be calculated by various methods. For a comprehensive 
ft distu^i(^ bf ftieflse methods, the interested reader is referred to articles by Ashley and 
II Ro^eii (rePf 1)^ LairfaM £tad Slhrfc (ref. 2); Ashley, Widnall, ahd Landahl (ref. 3); and to ^ 
• fehe%nore«rebeHt work of Mortnd (fef. 4)! The particular method of the’prbsetit stilds^is’h ’ 
cwnbinatioiv of two previous meMiods: (1) the 'classicrf ptessure-noftnSlwash meftidS 
I osiginating with Kiissner (ref. 5) amd developed further by Watkinrf, Rtinyan, Wo61ston, ' 
, Cunningham, Hsu, Rowe, and others (refs. 6, 7, 8, ’9, and H)); aridt2^he splined Ideal ’ 

■ basis, functions recently developed by Mm:cer, Wefeeft, and*otheik (Ifefs. M. and 12) for th# I 

■ case of steady fiowt 'Riis ftnmulatien using locaft basis femsticmsdn^ofve^tlfo An^dd/ I 
^Kiissne% integral eqnatiow oombiaes^ the smoothness and aiccaracy oA dtstributfon i 

me^|h(^dsj tly §igu )lyi(y und yergpt^i^ ^ panelMmetbcd^. ghe gregeuft ns - | 
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^ applied tp the solutioq, of subsonic, con^ressible flow problems about oscillating: planar ’ 
j , sur^c^s jt^itjb ^i^rqls^l^q^g^ t4>o^9tl^i])^t^a4 n)ptl|tp<^ ^Oloyed apply ta » largar 
^ clas^ q^8Qju|^oy^u^ y yuj^i^e^npyplaiy ^u|fa^s^et(g. )4 


I 

t 


i 

s 


Let S denote an aerod 3 mamical lifting surface contained in, say, the xy-plane and 
immersed in em ideal fluid of density P and sonic speed a. A rectangular cartesian 
coordinate system (oxyz) is fixed with respect to the referential surface configuration. S 
is assumed to translate with respect to the fluid at uniform subsonic velocity V in the 
negative x-direction and to perform simple harmonic oscillations in the direction normal i 
to S and at a circular frequency of cj radians per second. Then, under the usual ; 
assumptions, the method of acceleration potential (ref. 13) may be used in conjunction | 
with Green’s identities to obtain the well-known Kiissner integral equation, which j 
states that the normal components of velocity of the fluid and of the body are equal at | 
the surface. Mathematically .TCussner’s equation* • * ♦ I 


w(x,y) 



K(x - 5,y - n) p(5»n) dA 


(1) 


1 

I 

! 


(5 »n)eS 


is a linear transformation from pressure p to normalwash w, where the in^ici^d surface I 
integration is in the sense of Hadamard (ref. 14). . The kernel of the transformation 
, (ref. 6) is given by * ^ ^ 


K(x,y) = ^..p (-Y)(-*^[k,W) 

‘ H I. W) - H W- i] 


x-M r 


- 1 


/. 


THTf 




exp 




TdT 


VT 


+ f 




(2) 


"where *M j|s'V%sl^ number, r^=T^x^+?'®y* ,^ni ^^eTi j’^rid mo'3i^(FliSsel 

functions of JhCj first order aqd L i is the Struve fuimtmn (ref. 15); 
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V > 
> 


The solution to equation Ol) proceeds by expanding the unknown pressure in terms of a 1 
prescribed set {Pn } of basis functions j 


p(5,n) = *4PV* ^ b^p^(C,n) 


( 3 ) 


and determining the unknown {bn| by satisfying equation <1>; either^oUocationally or, ! 
in the complex least squares sense, on a set {(xni>3bi^of normalwash control points: 


w(x^,y^) 
m m 


Nw X 1 




(5»n)eS 


N X N , N ^ N 
w p w p 


Np X 1 


where and Np denote the number of normalwash points and pressure basis 
functions, respectively. 


f I 
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3.0 SYMBOLS 


I ■ V 


'**** •#* - r m SCa 

1^ SyiAbdls tis#d are listed below. To the right of the symbols are listed their definitio^, ^r 

^ the %ection, equation, or refexeace in which the 3 bai^ 4 eiined Of first appear.^ a 4 a 

L * ■ 

p Care has been taken to not use the same sjrmbol for different pieanings. There are ^few 
?■- exceptions, but it should be clear from the context what is meant. 

[a Any dimensionally consistent choice of units is permissible— dimensionaf or 

nendimensional. , , * • < 

■ • • - * a » «« • ■ ? ( 

t , A,A a a .- a • ^ preas of primary and secondary surfaces, respectively* * ' 

I “ ^ 0^ m m m > ,?ectioji A. 6 ^ m m m m m m 

[ AlCgjj^ ,AICjj^^ ,AICjj^ aerodynamic influence coefficient matrices of chordwise 
^ pitching moment, spanwise rolling moment, and vertical lifting 1 

, force, respectively, on panel 8 due to normalwash at point m, 

i equations (A-23) through (A-25) 

I 

Aj-b generalized aerodynamic force coefficient matrix, 

equation (A-26) 

i . ■ ■ ■ a a - "'«-'• • 

^ AR* * ^ • aspect ratio , aaftaainttM 

W * • ■ . V , r . ^ 

TtTo. 

* [a,bjja,,bj| ' * '* spanwise intervals of integration, equations and,(34) ^ ► r » 4 

. '* * 

■ bn * * ■' * » a » pressure basis function coefficients, equatiop (3) ^ ^ i 


f* Cl .^Lh •••**' lift coefficient for prigiary ^nd secondary surf^pe^^, ^spectiv^l^ ^ 

I • * > equations (A-12) and m ^ 

k* « 4 « 4 * • 

I • • • • derivatives of lift coefficient with respect to an^le of attack and 

1 . 9 -i ~ ^ angle jof, control deflection, respectively ^ ^ ^ 

? • ■ - ■ V- ' i *ytT 

i c.Cp ^ ^ chord length, flap chord length, fi^re 27 ’ ^ j 

c(i 7 ) , ^ panel chord length as a function of spanwise coordinate,- 1 

I ^ a IZ _ M 

'• I Cmiw ’*c3lfifci?hPn3&trix relating normalwash at point m dye to nth 


. _ pressure basis function, >equatioB (A- L8) 


^’9- 9 ‘ a • «HP • 

* • « 4 ^ 


cj(tj1,cj * spanViSe *^sectional lift -diaftrii)u|ion ^r^rjp^ and seconde^ * 

^ ‘ *sllrflcdfc, fespeetiuelj^ qtUgtlfn^(4-8iapdlA-.9) „ * * ■ 


vertical force on panel £ with dynamic pressure factored out, ^ 
equation (A-20) ' 


i 
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A 

G*(x,y,T?) 


.Wy ^ I • ' . * ^ j ■* *^*'''* 

stream.yiee integral, used fcg' calculating finite part of infinite ^ 
i integral, equatioiu(7) . ^ ^ * ' ' * * * ’ 

dervva4^v£ of G with respect 1to its thfi-d*co6rdindte “ ♦ ^ 

, doutbly ^ifferei^ced, Cauchj^integrable form of G(x,y,J?), 

^ equation (8) ’ . ^ 

G*(x,y,Tj) with logarithmic singularity removed/ figure 9 


• ^ , 


i-mM 


! > 

I e 

f « - 

I ^ . 

mJ 


I N(n) 

i . 

**^'r 


modified Bessel function of the first kind and first prder, | 
equation (2) . ^ ? 


modified Besseh function of the second kind and first order, 
equation (2) 

reduced frequency, cj£ /V 

Struve function of the first order, equation (2) 

■ characteristic length; also panel index, section A. 6. 7 * 


Mach number 


pitching and rolUng moments acting on panel £ , section A.6.7 


.j normalwash index, equation (4) ' 


4 

.M 


the set of indices of basis functions whose supports nontrivially 
intersect the panel with index n, equation ( A-6) 

the set of indices of panels contained in the support of the nth’ 
basis function, equation (A-5) 


Np,Nw 


number of basis functiona* £md number,., o£ nori|}alwash points, 
respectively, equation (4) * t ; ■ „ 




%Ns _ 

* f 


number of loading panels, section A.6.7 

■ i A ■* 

the numbers of rows and columns, respectively, of panels ' 
covm'ing the surface i* .- , ^ 

■ ■■ • , 'I 

^ ^ ^ d i 6 M » 

basis function index, equ£^io^ (3) » » » « j 


pressure basis function, equation (3) 

i4 %T ' 


loading panel £, section A.6.7 
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* ■ 


prtmary and secon4ary *surfac^ planfornas, respectively, *^1 
figure A-8 , , * * - ■» *‘1 


R, 


L V • • 

h ♦ 

U ^ ‘ ‘ 


support of the nth pressure^ basis' function, equation (A-2> » 


« * * - 

, .reflection of Rp, figure A-2, • 


5 ‘ 

t • • 

*v* 


# * 

, * 4 


aerodynanric surface ^ 
freestream velocity, equation (,1) 


• • 


• r 

• • * « 
. # ♦ » 

r • 


■* • 


.. » 9 *■ * 


m 

• • 


diagoneil matrix of least squares weighfiim fi^tors, m 
equation (A-2) *^~****' 


fluid velocity normal to aerodynamic surface 


t t 


r5i 


ft 0 • 

(x,y,z) - - 

(xm.ym) , r , » 

y ^cp * ^ . 

♦ 


characteristics function for support of the nth basis ffinclEion, 
with a syrmmetry condition, equation (A-4) 

* ^rearnvme basiaifunction, eq^tion*^) ^ ^ 

streamwise, spanwise, and vertical coordinates . * 

. it ^ ‘ 

leading-edge coordinates, figure A-6 » 

fi normalwash points, equation (4) * • ■ 


m ■ * “ 

,5 streamwise center of pressure, referred to local chord. 


« • 


k 

nyo'- 

► - 

k Y,; 


„ „ -r • - • , * * , 

^0’^’Y2 ^ spanwise coordinates in fairing, figure A-7 a • * 

* ' ‘ . a • ^ ’ 

0 ^ * spanwise basis function, equation (6) • • • 


POreek Symbfel? 


V m 


S » r ** 


« * 


a . f * 

t * * ' 

L 4 • • * * 


» - 

angle of attack, radians 

* 

1 "- M*, equation (2) 
pressure differences coefficienf, figure 32, 


r ACp pressure differences coefficienf, figure 32, a • • ’ * ^ 

► 

“ fl£^ aefTectio|j ig radians, secfSon 3.4^-^H*®** 

^x"-k*»»***’*****‘ *•' ' •'la*****'* 

® * Vafring wicfth^a^o^flgune A-7^ 

> ■ 9 ^ ^ ^Ftieal%patial c^dmate, seotioB A.6.7 « 1 A V * 



; 0 

. m f 

« •» 

.n 

s. 9 ^ 

«• 

« # « 
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horizontal spatial coordinate, section A. 6. 7 


spanwise center of pressure for primary and secondary surface, 
res|)ectively, equation^ (A-16) and(A-17) 


left and right planform extremities, figure A-8 


left and right control surface edges, figure A-8 


streamwise Multhopp angle, equation (35) 


sweepback angle, section 3.5 


rqw index, of the panel cover, .equation (5) 


Mesa function, equation (2?) 


column index of the panel cover, equation (5) 


streamwise coordinate, equation ,(1) 


local panel coordinate, zero at forward edge, unity at, aft edge, 
equation (21) 


three-dimensional spatial coordinate, equation (10) 


streamwise center of pressure for primary and secondar 
surfaces, respectively, equations (A-14) and (A7l5) • * * 

f ^ ' 

forward, aft, and midedges, respectively, of a surfkce ^ 


forward and aft edges, respectively, of the' secondary"* surface' 
figure A-8 ' . • - ? » v i 


density 


- dummy variable of integration, equation (4) • • i. • • • • 

^ 

* * * ' indefinite integred of the mesa function, equation ^8) • • • • 

•••***■*' m m m. m m. 

It 9 • * indefinite moment of the mesa ^unction^ equation (09k • • • • 

0 s 0 • * * * ’ '^ ' 

^00999* - freq'uency of structural oscjllati^n„ r^d^ans pqj: 

^ ft. • • • * * ^ equation (8) ^ i 

9 • 


0(»<y)t 


f 4 ♦ a^oByhamid panel, figure 1 

9 s 9 * ^ f y' * i 
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4.0 DES(^PTfoN OFPrEsENTMETHOdT 


* 


B. • 

• - 


• * « ‘ • 

- * •• 

4.1 LOCAL BASIS FUNCTION METHOD 




« 


* * 


4 & 

• • 


; « 

- • 

^ # 

* m 

# • 


A major difficulty with the kernel function methods has been in the selection of efficient 
and general pressure basis functions. Early solutions were as suggested by Kiissner, in 
which the basis functions were products, p„ = X„Yjj of chordwise functions selected 
according to two-dimensional theory and spanwise functions Yn selected according to 
lifting line theory. However, as analysis planforms became more realistic, especially 
with the incorporation of controls, the need for a more rational approach to the selection 
of basis functions became increasingly apparent. A major improvement was provided by 
Landahl (ref. 16), using the method of asymptotic expansions to exhibit the precise form 
of singularities in pressure that must be present along planform and control boundaries 
for the necessary jump discontinuities in boundary conditions there to be satisfied. An 
extension to Landahl’s mathematical analysis in the case of swept hinge lines has been 
pointed out very recently by Rowe (ref. 10) and has resulted in improved accuracy and 
improved agreement with experiment. 

The pressure basis functions of Kussner (ref. 5); Watkins, Runyan, Woolston, and 
Cunningham (refs. 6 and 7); Hsu (ref. 8); Rowe (refs. 9 and 10); and others may be 
termed global in the sense that the basis functions are defined over the entire planform. 
An inherent difficulty with the use of global basis functions is that the presence of local 
planform features such as controls, etc., affect the formulation of each basis function 
everywhere on the planform. An improvement may be offered by the use of local basis 
functions, which are defined over relatively small subsets of the planform and are zero 
elsewhere. 
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wner* is a|lojal, chordwiee>!CQ#rdin*te^ where is a cwtipu^u^ ] 

linear spline the local diordwia# coocdiaate, %md whA-e each^Y^ is a coni 
^ifterential&le^ piecewise quadratic ^plina in the'spbn^i^ codt'dfliate. Sinc^ 
fudctionS and panels have the same, index set, the number of basi§ fSnctions ei 


wise 

msly 

bapis 


I j’Q|e^pport o^a^n^tian is the closure of th^ sef on which the function is nonzero. 


The transformation from single' to doiable subscripts is given 


n = V + (p - DNgj 1 < u < N^, 1 < p < Ng ( 5 ) 
^It^^o^vqnifn%tqpiae liielBdhi#sA)^$t ^qfti^^ce' fOT^^e^ijj^ ^c|^l ^asis f ' 

^as^s^ajal^e^r^iicta^ v . - • ^ 1.1 i . # » • f - f f MULl 

P„(^,n) = X, (x)Y^(n); l < p < N„, 1 < v < N„ ( 6 ) 


* • 

) rows (a 


♦ 4.A Crib^DWISE FUNCTIONS t «■ » a • ^ 

>••••»"<• ,| 

Por*a %i^n'* value of y; the support* of is (1) tow hi of panels 4f p = 1, «nii •) rows j 

^ p-h aAcbp of panels if p = 1 * 1. The chordwise functions are depicted in^figpre 2. Xj is ^ | 
V defineck over* the first row^of panels and has unit value at ttm panel lespUpgypdge .pnd 
Vzm'o value ah the panel trailing edge,.Por p > 1 each Xp is defined over t 9 ro^a{iej^ row|, ^ 

• has zero wlue at botii the leading edge, of the front panel and the,trgili{ig edge^of^the 

L ;,afr pan«l, «n 4 has unit value at that interseijtion of the frqnt,pqd rear p^els^ « m ^ 


ROW 2 


ROW 3 


ROW 1 


Figure 2.— Chordwise Functions 


1 panels 


^ number of panels and, therefore, .thp nuipbar of norm^wash control^o|nte muat unlK :| 
I ^ |u^e<|tijp ^ 


■| 
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r^.r2^^VNWISE FUNClibNs' 


The support of is columns 1 and 2 of panels*if y = l\ columns 'P'- 1,1? and i»+ l**if ’ 
1< i^< Ng; and columns Ng- 1 and Ng if i>=N . 'Kie’fepanwise f&nctiftis’ai^'depicted in* 
figures: ^ ■ J n ^ 


Y,< SYMMETRIC) 


Figure 3.— Span wise Functions 


Eacifi Yy has unit amplitude and is continuously dffierentiable across panel bound^ies. 
If Ki/< Ng‘ Yj, has zero value and at the left and right extremities of its sapport.\If 
Y„ has zero value hnd slope on the left and zero value on the right*If*» 
has zero vahie and slope on the right and zero slope on d;he left^i&that edge, is. 
symmetrically reflectove, but zero value instead if -that sedge is aatiayiumetiic^ljji 
reflective or free. (This cSsesis««ot illustrated ia figure 3.)o a r 

■ aaaaa* 

4.1^ OBSERVATIONS ON LOCAL BASIS FUNCTIONS , *••**»» 

a-as'* *- t n * 

The above compl^ely defines the pressure basis fvmctions^sed jn the p|e^n^ I 

typicakpaessuos basis function is shown in figure ^^r^ ^^ryatipq|| ^^b|^ i 


^ in r|^^d^ this spt ef ^n^ti<^. 


Figure 4.— Local Basis Function 
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f . ^ I? i *a basis becai(se4t is UmearFy fe^pen^i^ _ _ 

• ^ 

2. • If incomplete m the Jiipit^as the maximum' panel dimenaioo appeoaches HerS. • * 

3. - It nfey'be extendi topologically to nonplanar surfaces. _^,##*^*'***** 

• ^ < 
4n A Finite '’suins, used in prentice, , provide polygonal a^prox^aticpis, cbondwiso'aAd^ 
^ » qua4ratic*appr6ximations spanwise. ^ ~ * ' ' ' ^ *,.«*#*^ 

• ^ 

^ g|F«irtheT %(Minons sucl^ ^ ^n§oi:f>oratto4 of^tfape factors #natke<i 

H uLignkaM-lloMfe Sin§ufai^i^ ^qg))%m§da imaktAi^tforw^d mapi^r^ • • # * ^ 

g E^eRA'g(^^^q|I\Dk£«Ifil(3Ue^l^ ^ 

The coefHcient matrix of equation (4) is evaluated using Gaussian quadrature as in 
reference 17. In this operation the kernel function, equation (2), is separated into 
nonsingular and singular terms — equation (9) of reference 17. There are three singular 
terms identified as dipole, square root, and logarithmic singularity terms. The 
nonsingular terms as well as the square root and logarithmic singularity terms are 
readily evaluated. The dipole singularity term, however, requires special attention. 

• • • • r , 



b 



Gft(x,y,n)dn 




+ G,^(x.y,y)log| y ^ 


• wkdie *llifl lLl£cyp^^^|da!n^ai p»rHa^ ^lf(^efti?tim^v^e|^ ^ ksMtlft ® 

Igliate^rA dt tffe*8ipole singul^tli tvios th^i€sArf, &(T where maSV*****^ 
— — -■ -- - - — ^ T ^ 


s.'-.y-”) = y-h; 


f) - G(x,y,n: 
y - n 
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p* The coefTlci^nt*iMtnx^^quaT?oir^f^ i| b^si^ %ivibiwi4t •« 

I" index and the normalwasfi con^rof point index. %o^e^e|^ when ]pc^ ^oi^tLona a<« » » 
r* employ^dt it is preferable to index the cal 9 ul^tipn^ pa{)el^ i^tt|ei^tl)p% basis 

functicAis'’ to avoid nnhecessary repetilions ‘ in calculating th^ A{plije^o^ tj^e. kernel • «4 
function, it is consequently necessary to determine exactly which basi^ functions a^ ^ . j 
: '■ supported by each given panel. The required transformations between indices ipay Ijp . 

I * obtained in a straightforward fashion from equation (5). Also, it has been foun^ ^ • 9 i 
cflieficient matrix inversions using local basis func^ons have been well condition^. _ « » ^ 
. A - -a 4^ m A A A K A- A A. wk S ^ 
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«.e eOMPTJ'FER • • 


^ tCjm^utera-eflulte *r^res^tfed%#sil e^vS^^impfe slea^^si 

^ «ii|l A^ieh «ttic interiiiedKaAe ^iStHfenfetfcsrf A;cAr^cy* evafuale^ a 
k <alA> *th6 fhtfriRedlafe Vi^tlieinatlc^ 

I ^^wl^ated;•(3^a4DOBsteally f>r«bltn#ir(*wfti(% A>nlfei^efkc# i^sAdfed*ai% &ii 


^ «ii|l A^ieh «ttic anteriiiediaAe A:cAr Jb; 


accurac; 


le results are ^ ^ 

comparetd wMjh .lanother^ c®de;* (4) a noristead^ dorilrd! SurfacI prot>le’m*ln^ wiii^ H • ' 
rejuljts ^^ccgnpared with anothar «ode;45)*a ferte^of^^nfindtrft p’lajfofmS consisting 
ofjSkewed p^^leIo^ams]| and* (6^ the freqij^nKyfreepen^ fif % la^ef Abfoifcc^eF * * 
^r^spoi;|t to ^s^lqpdiPg. I)rofaieins}(lhand (2>lahB bbsfc <Ae&oht)^rdbliniS iA A^i^ ^ ^ 
^upi^ic^l (espH^ ctieoked by aeparate<calcuiBtloift, Wh4rdhs^r9bt9m8 ( 9 ) 9) 1^) 8rF ^ ^ 
pri^blfm^ dps^|^qpd,it<\|V^i(j[pta the^pDo^iaia by voapariaoi* ^thf the %-esuIls 9>f%tfle^ 
methodSj^ I}ie„cqfni^t$r £0<|p vs^ lor^the present otudyais^ numdbicAil 9e^dSctiP 9 ' 

54LiA‘f ' 

^ ^ m ^ ■ m * - -ibk-f, 

5.1 CHECKOUT PROBLEM ONE 

» ^ • * • " » » • *■ ' ' — --r — ■” 

P)pi)pa§) t^ie iimt thackaufrpaoUem, €Pl, fb ik> (Irc^ift 9n%cAiAtF Ao*cAbffn* ^ 

/>/\mr\Q v*i o/\n ark^aci LJ n Ml vJ-ilm->-<i <4K/Mb\i4lVktftr>9/1ac/^KTi^(rMr 


Qixvxk.vuv ysuutciU) ro vu pxvrvxvro cxxx ovv>uxau^ x^xuo^vxtxuxax 

^nj|p^i£|p^oi||tl^ ^iqplagt 4poasible<f>rabknm 'Bke^plstoftrift dh(^i^bAd9dftmpfficin 


Simplest Possible Case 

k = 0 2 panels, run as an asymmetrical problem 


mid panel downwash points 


n 


DWCP 1 

DWCP 2 

+ 

+ 

Panel #1 

Panel #2 


uniform angle of attack 


Figure 5.-Planform for Checkout Problem One 
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' 9 ■ m ^ 9 wmm 


! For'this problem, thp coefficient matrix^r 4 

^ J- , 


1 


1 1 


't 1 


I a 


'a a 


l::l ■ l::l 


l-n 

is doubly symmetric. i,e., 1 

j * r* ! i n * .1 


C j j - C j j & c , , = c 


1 a ''ll 


and for*unifonnMoVriwash\ 0(^ «2 »a) i^f^l(j|ys,that 


b, = b, = 


8 ir a 

C — 

II ^1 a 


I Thte «u<|Wi’^te^fc^l®%A 'fh^ ar^iven by . 

L > ma^ III T 1t " r> t M 

c =-t’^ 

mn J 
o 

wherQ;-"! 

- MaJ t 

Sn'-m-yn,-"' = ’'n"»/'l ‘ 5>[l 


(y^ - n)* 


"m - ^ 


(9) 


( 10 ) 


(nj 


( 12 ) 


Vcx^-O* + 


(13) 


and 
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The first point of comparison between CPI and RRX is It may be shown^ 

that * " i 

G(x,y,n) = + (1- |x) + B*(y-n)* - ^(1-x) y< l-xj* + B* (y-n)* 


3^(y-n)^ log V(l-x)^ ■»• B^(v-n)^ - (1-x) 1 

* X + -Jx* + * J 


06) 


yx* + B*(y-n) 

^ where, for simplicity, the? subscripts have been omitted.*It is necessary to define G(x,y,y) 
as a limit; then I ^ ■ 


G(x,y,y) = lim G(x,y,n) = (2x-x*) Y(n) 

n ■►y 


(17) 


The graph of G based on equations £md (17) is s^own in figpre 6>A 'snfall "ad^oc ~ J 

4. i?1j 1 «L_ ^-1 '' ? T 


r 

1 digital computer program called ECPsJ vnas^-used Ito calcula|n qiPmerieaP<vdlues forlall 
' the closed-form resirftrf*^reSented^n thi|^sect9kn.'Thes^vaIues,^ej(e tben^ cSPinffered with 
the myn^ical values^btkined’*by RRX., The «eillie#obtaine(f^y RRX for (i-wnirefoiliidln 
be accurate tq about seven and one-hajf decimal places with an average«deviatlon iH 
about the eighth place. ‘All spanwise integration statJoh's for rj were compared. A'^brief 
► excergtw^^^^ in table 1*. — _ v * 


Table 1.—G(x,y,T\) for x = 0.5, y = 0.5 Basis Function No. 1 


n 

G(x,y,n) exact 

G(x,y,n) RRX 

♦ ♦ ^ 

Error, u% 

.0267 04S5 S463 

.5008 8326 9740 (-1) 

.5008 8326 8670 (-1) 

-.21 

.4919 B863 3361 

.6982 5165 5022 

.6982 5169 2610 

♦ 3.95 

.5060 1136 6639 

.7072 6110 0380 

.7072 6113 3184 

♦ 4.64 

.9732 9544 4537 

.5030 2766 2341 

.5030 2766 1266 

-.21 

1.976 8560 5193 

.2227 9916 0683 (-3) 

.2227 9916 0793 (-3) 

• 

♦ .05 


*^C0M 


a^tations Filg^ Run^CPl.6 a^d CP1.7 
** ’ One/i% e^oj equals^e^art m 10® 
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G(x,y.n) = ^^(n) 


1 + (2-x) /x*+6*(y-n)* - 0-x)/(l-x)*+6*(y-n)* 


+ 8*(y-n)" Inn - (1-x) 

Vj^+8*(y-n)* + X 



Figure 6.—Chordwise Integral- Downwash Point One (CPI) 


r 

► 


’Tffe*inft^irarfd fh^s^riWile Integral de|«cted^il*figure 7; it suggests clearly Uiat 
th^^tdfral dod^ ffol^:dsl^^a# ordinal totegMlWeirtliep, iu ouder tu uiew the«fiect 
*of^VI6 trdhsIbrftisAion, th# graph tf tbs' singly idiflierenoedcfuiictaoik is presented, in 
♦fi g ' tu O 8.» It majt 0% thah thi» difference' funotioB diSferg %o« its ^ limit, ibe 

deri^tiv^'G,^(n9y,y) the oed* ef |y-t}|log|y^|neer« 7 zy«eq|l Itia^ ^ infinite ^loj|e 
will>be exh^iteddn*fl|yute d ah the point nnd, furtheu, thfd t^e^ohbly 
hinetian G* wilNpoBsasswanlogaaithmie ainguleritypthpr^ is 

* integwhln, as ahw^igiib figure 9. KemovaU-ot it^ Igg^ithmic d|pi^tM in 

rigurapJiO;^he hvietioct is completely nogpuig^Ey;; a^d }i^c^n^y^^q|ily i^^t^pated 
^bytstend'mi»hig|li~prtpi^n,Qaus8ign.m ^h^|^ ^ ^ ^ ^ 


A 


' I 

« I 
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G.(x.y.n) ■ G.„(«.y.y)]pi^ 

*'mn “/ ■*■ G,,,(x,y,y)log|^| + G(x,y,y) [b^+y!?] 



Figure 9.— Doubly Differenced Spanwise Integrand (CPI) 
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'n 


I 

f 


G(x,y,n) 

<y“5* 


dn 


m 

■! 


G( x , y.n ) 

(y-n)» 


dn 


( 18 ) 


^ We the jn1tog/9tion ^lia^part in^lvTng^the 

analytic 4»tegrand. From>«th^ standpoint irf rfRmeripal analysis, one. may use differ^t< 

» numbers, say NG”’bf Gaussian quadrature points^ aq^ different numbers; say ofi 
„ subinteryal^ These respective "choTc^ ace indbpend^t.-and com^r^ive c^fcuKatmna^ 
^ were mad^ far NG = 2j4,846*<an'fr NS = 1,2,4, rasults 9re p^seotec^^ bflow ^ 
table 2*'^e^6^iht singleiasubinter^l pontfaiiUtion^conesiionSt^ liiG ^d*^na 
„NS= 1 appe^s<(^^num^or4bi»caSe. It jis clear tlial^onjtprfende is ^abipde (lC*may 
|^ba<ncifSatl^t cipr^re^id cqgyeageltce m^ 1)^ obtafted b^ufipg%he’*d^l4K difTdftnrac^ 

rinte»rariaG*.4.^^2^ ^ 

^ 


Table 2.—Legendre-Gaussian Quadrature of \GU ,y,v)/(y-T})^ dr] 


NG 

NS > 1 

NS > 2 

NS « 6 

NS « 9 

2 

.2137 7t2S S«36 

.3166 9639 1676 

.3219 9629 0633 

.3226 1336 6362 

It 

.3220 tSOe 666>l 

.3226 6332 6213 

.3226 6611 7621 

.3226 6622 9660 

• 

.322S S«20 3769 

.3226 6623 0226 

.3226 6623 0263 

.3226 6623 0263 

11 

.3226 6623 0266 

.3226 6623 0266 

.3226 6623 0266 

.3226 6623 0266 




M 


1 1 

f ^(yini? ^ ^ 'll ^ G(x,y,y) + A 

0 o 

♦ G,j.^(x,y,y)log 


( 19 ) 


* ComputationrFneTRun m 
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IS for purposes*of numBrical inte^ation, a o£-integrating G« for • This is 

conveniently ‘done irf tiWs iJhrfe:iie (^0,yj and J siAply uses standard 

Legendre-Gaussian quadrature, the numerical results may be expected to be poor 
because^of ,the4ogarHl(fmic singular)^ remq^nii^g^ G«, thaU-is displayed id figure 9 
Mg iSgthft cqgejeas «h(^m in thfole S.^ow^eigj^ r^jp ^e,^mt«lysis peecedifl|r fltbntirichl 
integration, , the logarithmid singiilarTly^ ^fparated as it is in figure 10, then 
Legendre-Gaussian quadratue may be used on G« Euid logarithmic-Gaussian quadrature 
may be used qp the 4pgaritlmitelly singular part. This i^lSA shovekinKtaWedi UNnay 
be seen that the eight-point Legendre-logarithqiic form of Gaussian quadrature 
converged to almost 12 decimal places.** 


Table 3. — Gaussian Quadrature ofj G* (x,y,r\) dr} 



Legendre Only 

Legendre- logarithmic 

2.512 6234 7443 

2.770 2077 8297 

2.582 1430 0494 

2.770 0346 6808 

2.603 6269 2847 

2.770 0355 4307 

2.609 6812 3399 

2.770 0355 4308 





1 The cumlilative convergence of the numerical integration proc^ures for C mn is shown j 
in table 4. It is seen that the convergence is quite strong. The level of convergence is 
; between 11 and 12 decimal places. _ _ ^ ^ .. JBe-.-aJ 
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j By comparison, the results of RRX are shown in table 5 to eight decimal^lacesTlt is 
. seen that the maximum error is'0.10% in the smaller, o^diagonal term and 0.04% in 

I the larger and dominant diagonal terms. Because the integration procedures^n RRX are 

i rather extensively modified versions of those used in the RHO-III code, and^th'e latter is 
I optimized numerically* for g^abakbaais functions rather tha^localT basis'Yuhctions, the 
^ agreemen^sttowif in table 5*is ^ifficient for ovif p ^^rp^es^ ^ % “* 

Table 5.— Comparison ofC^^ From RRX With Values 
From Table 4, Which are Essentially Exact 


Source 


Exact 




5.259 9813 


*"12 ~ ^21 

-.3842 0736 


5.259 6758 


-.3842 6833 


Although a two-panel solution is extremely ^crude, it has Feen most interesting to 
compare these results with those of other codes. On ^^e Basitof'the above results^ it is ' 
easy to determine the, basis .coefficients and then to compute Ci (which is a gross., 
measure of the ability to predict pressure accurately). This was Sme, and calculations 
were also performed (1) with the RHO-III code using 48 comparable terms (actually 24 
terms per semispan), s^d (2) with the steady-state method Sf F.*T. Johnson (ref. Ihl-taThSk 
comparisons are shown m tabled. The agreement i#rdinai%a1ile. • & 

^ ^ ^ ^ A , 


Table 6.— Comparison of C / From Different Codes 

‘-a 


Method, of Ref. 18 
RH0-III 


100 terms* 
48 terms* 


52.46 

2.474 


2 term * 


.2.577 
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^ TheToc^ panef co^dma^s^e then giveii,by*. 


(VC: C, i 5 i C,)<r = ^ ~i i ^ ^ * 

_ ? - 1 “ ■rl*il 


( 21 ) 


and the funetioifs are given by 


i 

i 


(Xj = l-I) & (CjiCie,— X, = X) * (5,5515, -i-X, = 1 -X). (22) 


The stream wise functions, X|i, are 'depicted generic^lly in figur e 2. For this particular 
pr oblem^Jh^panw^ functoon sf^^^are depicted by figure 12. 

(Oin<i— Y, = i-^*)«i(i<n<2— Y, =^(2-n)*) (23) 

(osnil — Y, = ^n*) & (I5ni2 — Yj= •j^(-io+20n-7n* ) ) 

4 (2 i n 1 2y — Y, = YJ(5-2n)*) (24) 

(li n i2— .Y, =|(n-l)*)4(2< n l2i— Y, = l- [l-j(5-2n)*]) 

(25) 


pointr^spe^tivaJi[, a^jg^frc^^en^ omittQ^|fe^w; howe^r, it shoura ^^ear^ 

««tUa 4- ia wkAAnf ^ 


|! doam^Msh 

from tl)£ coMtext what^s ntgant. 

^ -«i» ^ — 


G(x,y,n) = Y(n) j x(T) p(x-c ,y-n)dC , 

(C :p( 5 ,n)^ 0 ) 


(26) 


p *TH r^.p\ (r| satisfying^ th^|also satisfies O . - 
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♦(x-e,y-n) = J u(x-c,y- 


n)dC 


^-C»y-n) = j (x-Ou(x-c,y-n)dC, 


T' and* one finis 


♦ (x-C,y-n) = - (x-O - ^x-e>* + 6*(y-n)* , 


♦ (x-e,y-n) = - *s(x- 5 )* - Js(x -0 V(x-0* + B*(y-n)* 
+ ’s6*(y-n)*log [V(x-5)* +6*(y-n)* + (x-t)] 


Indeterminacies iS* 0 must be handled as limi^. In additio^ it is eSSential to ^ 

recognize that alohg tj = y, the behavjor of certain singularities* ^penc]^ bnljp>u^on the ^ 
sign of X One finds • ^ ^ ^ 3 

G(x,y,n) = {cej(n)-x][((i(x-C 2 (n) >y-n)" 0(x-^; ^ (n) ,y-n) ] 

(Vn:i<n<3) + [0( x-5 ^ (n ) ,y-n ) - <l» (x-5 j ( n) ,y-n ) ] } .02) 

G(x,y,n) = I Cc , (n) - x][0(x-c j (n) ,y-n) - ♦(x -^2 (n) ,y-n) ] 

(Vn:4sn<6) + [0( x-^ ^ ( n ) ,y-n ) - 0 ( x-c ^ ( n ) »y-n ) ] 

+ [x-c j (n) ][0(x-5 j (n) >y-n) -0(x-5 j (p) ,y-p)] 


- C0(x-Cj (n) .y-n) - 0(x-e ^ (n) ,y-n) 3 j ,03) 
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«(i}) is the locat^wel chor 

► -m 

L ^ the abovet'Closed^orm expressions, the*“^lues' of tJ(x,y,fj), G(x,y,yJT'^Sit 

^ were computed’*' with an ad hoc digital computer *prograta calTe<r’E( 

yJ^pfy)utations were made at each downwash point and at each of the many spanwise* 
^integration stations used by RRX. Several thousan<^ numbers were compared;*'" th€ 

; results of RRX for the CP2 were as acc.irata as in OP1 . > 






^i^ls. The^ 


dr this problem, Cmn wi iT d^end^ upon integration over as many 
dwnwash point will belong to, at most, one panel, and the Hadamard integral must her* 
ferformed for that particular panel. All other panels support analytic integrands. Thus, 


mn 


E 

(a<y<b) 


[/.' 




dn]* / 


G*(x,y,n) dn 


where the' summatiomiis performed over th^ spanwise coordinate projection of those 
panels in which the<«.dow.nwash point is not contained, and the last (Jjre^ terms 
correspond to the one in which it is. The analytic'int'dgrands on the intervals integrated 
by ordinary Legendre-Gaussian quadrature and the logarithmically singular integrands 
on I a«,b,^ I are integrated by logarithmic-Gaussian quadrature. On the basis of these 

*caiculations,ia;* the • . ■ ' - ^ ^ - 

■^Values from«if{Rj^8 


L - 


he ^ e^coefficient matrix '•for CP2 4s shown, t^^u^eq^al places, 
^ac^^o^* i^^arSh^l^s'helow thesexacjt M^ljue^ ^ ^ ^ 


IT" 

1. 

h 


Vi. "Yalues^f = j^must be worked out anj^lyti^lly^^m the aBbve^fexpr&ssions 

due to indeterminacies. 'Limit caleulatiansvan Urn computer using t!CP2 we?e nRalljnill 
condi^one^ ^ ^ ^ ^ ^ 

•*RRX Calculations File, Runs CP2. 3 through CP2. 9. ^ 

***RRX Calculations»File, Riu^s CP2.10 through CP2.1^. 'Ae r'Jfeulfe- converged toil decimal 
places. ^ ^ ^ ^ ^ 
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2.41 40 S 

.729 411 

-.034 795 

- 1.76 993 

.066 543 

--.044 384 



( 2.41 419 ) 

(.729 438 ) 

(-.034 795 ) 

(- 1.76 961 ) 

(.066 607 

(-.044 384 ) 



- 1.23 159 

3.58 102 

1.12 170 

- 2.26 456 

- 1.58 112 

.186 115 



(- 1.23 161 ) 

( 3.58 089 ) 

( 1.12 167 ) 

(- 2.26 454 ) 

(- 1.58 098 ) 

(.186 145 ) 



-.465 187 

- 1.86 858 

3.41 598 

-.822 622 

- 2.80 607 

- 1.63 181 



(-.465 187 ) 

(- 1.86 860 ) 

( 3.41 585 ) 

(-.822 621 ) 

(- 2.80 606 ) 

(- 1.63 169 ) 


[c ] • 








an** 

2.23 192 

-.679 055 

-.119 257 

4.69 607 

.028 633 

-.139 504 



( 2.23 192 ) 

(-.679 013 ) 

(-.119 219 ) 

( 4.69 620 ) 

(.028 760 ) 

(-.139 504 ) 



- 1.13 757 

3.40 261 

-.108 931 

- 2.27 143 

6.95 942 

.964 842 



(- 1.13 757 ) 

( 3.40 261 ) 

(-.108 930 ) 

(- 2.27 145 ) 

( 6.95 929 ) 

(.964 812 ) 



-.483 378 

- 1.75 293 

4.63 633 

-.924 643 

- 3.49 212 

8.13 560 



(-.483 378 ) 

(- 1.75 293 ) 

( 4.63 633 ) 

(-.924 643 ) 

(- 3.49 214 ) 

( 8.13 547 ) 




IK 1 • 

19.68 

476 

8.25 933 

14.8471 

2.25 404 

2.36 621 - 1.74 502 J , 


(9.68 

387 ) 

( 8.25 935 ) 

( 14 . 8476 ) 

( 2.25 420 

( 2.36 619 )(- 1.74 532 ) 



L " 

3.74 

249 

(3.74 241), 

a 




cp 

.252 

220 

( .252 224), 

cp 

. 432 

977 

( .432 980). 




3.742 (RH0-III using 24 terms) 
3.784 (RAMROD-X using 6 terms) 
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irwxi W’** 


*TlS*ag<fe Aa^R*^j <ii> six-pa^el^ 

quite crude. Since^hfefe^are onlyil'w^lifeSBmr^ row^orpanllB.lthe^cfia 
distribution consists of exactly two straight lines. This representation is illustrated in 
figure 13. Although no comparison of pressure was made with the RHO-III code, the 
sectional lift and center of pressure were compared and the results are shown in figure 
14. The agreement for sectional distributions is quite good, but not as good as that of 
the lift coefficient. In general, integrated coefficients agree most closely and are the 
least severe comparison, sectional distributions tend to agree less closely and are a more 
severp compa rison^ j.qd prea|stu»e compain^lns a^e' the^most sevgjg^nf all> 



Figure 13.— Pressure Distribution From a Small Number of Basis Functions 
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Figure 14.— Sectional Distributions From a Small Number of Basis Functions 


31 


_l 8-1/Z X 11 INCH CROPS 





^ ^ 5^ YAlJQAyiQNJ»^05LEM ONE 

V* f 

I 'Thfe c^ 

I f>r#>lfcaA)j^ (] 
p eu]^ort4s^e^reaB^,*i|ll 

E HeopgaizSd %-eliabilily 1 
dtw4>vlis6 nioife ic 

> this particular study, 

g gHO-UI code.^his co<^ wa^ adopted ip f968 by ^GARD as th^r 'prjpaarj^ standard for’ 0^ 
J ]^^ady aerodyn'aniics. « ■ ^ ' * 

Cal^ulptjjpns far VPl are arrapgdd ^to fju/ s^es: the GP2 series, A |erie^ rf s^ries.i 

CP 2 ^^i^\ ,T1^ first nonsteady runhnjde wj^tf RIlX" was the same as CP^ except with 
tfie reduced fr^ucncy set equal to unity. W^i^as f(5- the stqadfy -state case (k = 0),* tl^^p ' 
error in v^s4%; fti the imsteady caae (k = 1), the error was ^2%.*- 'Hiegiextent of 
tlfis*'chanee ^s indicative of J^he difficulties enqpuntqred as one passes ^rdin ^teady to 0 
ufisfeady flow. » ^ * ’ . i i* ^ — i 

• ‘ • ' . - ' . # 
T%e,*ll^st*ques)tiq(i addressed is what afre the ^provemNits that are^ ^i^e<^ by * 
dlr&sin^ tlie nu^bpr of panel^. It woqJd*be^natural'to expect a greater imp^o^enfei^ 
frthilSnSreasmg^th^ qpmber of choedwise panels than from an eqtfiv^eht increase intthe 
feiiHnber* of spanwise^pepiels.' Thei^fOTe, the**s^nd ^run used four chydiyiag ^anefs 

linitead'oflinly ^wo. 'llii&. caused the eiTof tq^ drop from 22?%) to 7%. ^ f f ^ 

^ ^ f 0 , » y>--^ 

^he aesUltd of tfie lift coefficient calculations are depicted in figure' 15 using the complex 
^latievrdpresentation. The pomparisoi^for the stea^-st^^ calculations is sh<|^n also, ^ 
y^F#^ari>g talong* the real axis»>'nie peypeka^ajigePient^ ^e indicated as ^^njijcnhey of 0' 
jj^hardp'ise ^afielS fimes thg number pf 4span\yisi panels. Thus,.>for examjje, 

Bn(Jic%te% faui^chorffwfse j)apels and three spanwise papels. This qame cepre^ntaffbn isf j 
^s^ fo% t|(e <IH0-Ilf calcmlatitms, but im ttfat^cefee^ 6Cx4S (Jenbtes lix, clfoi^wis^ 

I pressure, funeftons by^our spanwise pre^uae ftii^tfons. For' the sixipanel c^cdlatiohs^ ^ 
^he s^ticppl distributions for^ li|t coeiffimef t Cnd efen^r^tof' pigpslui^ are<tc^peredidn ^ 
£?igur^l^ the^ 2tt term resuks^on^ RHO-HI. ^he •agyedhient is ^eji to 0^ ^ 

r g^uantita(iv^l 3 ^ %in and*qdalita^ivel 3 p *ri^ck Additional compnrisoms, using fciye^ ^ 
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Lift Coefficient, Out of Phase 
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Lift Coefficient, In Phase 


Figure 1 5.— Nonsteady Lift Convergence— Validation Problem One (CP2 Series) 
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Figure 16.— Comparison of Nonsteady Sectiona 
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oalled'Miilthqpp e<|uian|pila^ partition (ref. 19^ in ^Htch Ae cfford^se> dfndp ^ 


^ _ 1 - cos6 


0 < 6 < TT 


(35) 


n = cos^ 


0 S ♦ < ir/2 


(36) 


— » 1 1 ^ 1 ^ ^ f ' ^ i p I wy ■ T ' . T*"* " " " '• ' ' • 

w{)iere^qaa^ an^Iar yicrements in d^nd 0 are%sed.*The'feffec**i8*to igcreaip th^ 

paq^l deqpity in ^e're^idhs of^arge^t prSssuitf gmdients; i.e., along the l^adiqg anc^^ 



ri^ th^iumBer orspan^iae panels pwducod ap erryr of gl.59ybut ^itH ^ 

more accurate 'p^ase angfe, ad^mayiHje seen, by cpferripg tqi,ngvu;e 18., On the^otKe^ 
hand, Soublfng iii!&tead'the number of chosdwis» pqp^s {educ^ th|^erro|^to 5^%. is J 
apparAt^that*in*the ^jreselSt method," it* is paefere^le tg. h^vp ^ larger *n^fe^o^ J 
chofdwis? pafel# €han spainwiso panels. JFi^^lly^ dousing Jjoth^ the number^o^ • 
1 ^' ^oftwis? and*spai^ise ‘panels seduced th^ e^or qven ^rth^ to^.5%>. I^niaynje J 
g (ftseftVed4hat tbe 40x 4Sit:ase ihasneai^y th^sam^magpitudp as t^e 4C^ ?S 'fese, Cut m 
^ tJftat rfle 4(T x 4ST case* has a more «ccur«te, pha^ angler Jmgoifant criterlon^irf' ^ 
^ i#nh1elhs siAh a-flutter, .and ia an affect-that .the comnlex error ma^nitutfe cr4erioif • 
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Figure 7 7.— Equiangular Panelization—VPI (A Series) 
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ti 


. <• 4C » A 

B dkrjesAl%el^«Vp«BaaoiM;liiB«eMea is isveetigalp 

RRX aode te faii«liaa1iBO]»d«fisitjiki]» terms of accurvjiiiav ^ 9 f^n||)^i^d 

with tihe A aerieassC^alrufetwns. 

mm*. 9 , 

'{h» paneliistviw v^»asai# ^ 4P -^ 14 ^ 

ho\ye^r, ^qfl ;i>q)U:^aj|ie§ ^^eguulitfaAt^a^ ii^V ^ # 

prapti^al ipij^jj^i |[^v^i\^^s^f||ii||p^i^j|iv^ 

th^e^i^D|pi\§r j>^e|jp;^i(^. youldJb«^^c^ t^at t^j^e ^^lytign js »Utiyew 


is use5: In particular, 2Cx2Shasan error in Cjl of ^ .4% , the 2C x 4S of 22 . 5% , the 4C x 2S 

** of6.9%, and4Cx4Sof9.2%. “ . - 

W * ' • • * » 

J^^e ?act tfia^^tl^ 4C X 4S ^nelSaflon is^^s ?c5irat^t^n'^th’l -ft x 2(§ ft afic^^oft. 

Once this anomaly is recognized, it is seen to be present in the 2Cx,4S (22.5%) and “ J 
* 2C X 2S (22.4%)j)air, though to a lesser extent. In retrdspect, it was also present in the * i 
' * A' series: coihp Se**2Cx^S (^l.S%fwth 2Cx"'2S (21.2%1 The' s^arfwiSe ^aRef Widths^ j 
"nftiAHfe AoFaife fcifth^jftro^di* in lRe%qfii<fi 8 tan<f cifee*h9n*!n4th#e(fu?ari§€ifcu4iEafe# f 
Since the collocation points are being placed at midpanel,* it is seen that, with ♦ 
equidistant spacing, the distance from the apex chord to4:he inner eoltecartien point»isB 4 
muA Aoftr*thaiAvith*th%> e^uihngvdhi*sp(tclKi^ Slta4h «heM 3 (pkxeai 4 hiAeate*i 4 |'vdb 4 t)^ J 
iff (ftwnwtisR tllaiP i^Ry^aH^ifadbiiesille^iM drrehoeuV, eellicAidk poiHt 4 |M 

dhA to %uih *% %i 4 lguiarttyr w<ilh produce* ArineisuR icsaltB. Vusther# tliie ^iisape 
opetator K SnA^r# t#b# dlhcMittnVouB for plakfimfe erith bveaksi lUiauiqsaneflt Ikat t|^'J 
av(||d*pl*ciAgtiu«n«fceKsary«rflStiictaonB on* tke 4oaatiiBn*of aeU^atiqtP 4 pqj|it 0 > ^(hgi a 
i|la|pfarni>naiati ba sciCDtlKdk andstht •maotifaing vmst^iiapraatiQa, ^ ipppl^l^ ‘the 

» kiilik'WrinaeldioandesieaawHall* a • • • ^ I 

s . . _ , 

calcufoti^s i^ to rpflpe ,|.h£ j^an^luation I 


f 

* 

i 


C Series: Thg ^ 1 ^ 0 ^ qf 4he. Ci. series of cal ^ ^ « 

fuAhIr and to exanune t;^exoa^aigmiale df tffe f r^Sre distribution. 44 ^ * 




« « 


The c<Jh>lfergence'^f^L ^ 8®®* ®* »glartee1h(^%g^r^l^ Using ^6,, and €4 

l| ilanfel#, errcns ^are §1^^ 3^.5®, «ind 4.i>%* r^s^tectively, The^ qpxt l^val of 
# cdmfiarlsdh Tb df ^e^se^tigpqj distributions, in figure 22. It^s^egp^hat 

. . J.\ ^ 4_iA. a fe J Al lx_ Xt l-STTr-V TTT _ ^ J _4 *_ ® ^ i.1 X. •_ xt 


II 

4 



cod^exfcept *ear th^r^t^n^the * | 
level Rf liomparison is^tl||it^f^he • I 

An >A^4 .T A t f* ■ 


, -_ shown in ilguree 23, 24^ and'25 for 4^1g, m 1 

i 4 >aoelB, respectrvefy,*for ?lie^ nud-semi^aji iherdr ft>r* a' pifelimina^ stgdjj, J^e, ^ 

I ^onvergmice n#a 3 f b§ cRnmdered ^i^ goo^. ijh* approxinfoti'onal ftxutations imposed 
I fcy cherdwiae WneaiTty^orthe funcj^iops^lsp may be seen: * * ^ ^ ^ ^ ^ 



requiremei 
■ *(Nw > Nb 
I •(dtherwfte' 


* ^and every. pancL and some nanels.may JiayeZmore then (me j:o| 


"A 

iomt. Tfie 

c^ntrol^oiht? effuaT oF exc^d^he n'hmber of panels 
hf cfti^ol^pt^ltf ^ ®Aid^eAc®y%ifttlft ftuSdAyfof’any panel • 


l^y^pan^ s^n^ 

isents are that the nuiAer of con^“' '■" — 

tlRt'_ ^ _ 

t>f«!oti>!i ^ * i • 

i£surf SLM A 
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Figure 19.— Equidistant Panelization—VPI (B Series) 
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5.4 VALIDATION PROBLEM TWO 


a The plffiform, 

hcoefQaiaBLU-veij^u^£^U(;^^equ^^, amd chora wis^ pressu're^Estributioff'are shi 
Jmn Sgucea- 2g, Th g ^ control sui^ce a '*40^^^fioT'd full-sparf'*flap. 

equiangular spacings from the leading '’^ge to the'hfH^ 
Kour eqq iangi^lar ^pacin ^ s^ f^p m the hinge line to the trailif^ edge, '“and 
fcquianpplar spacings across the semispan. The agreement “with Rowef'^s^aeen to 
l^^lept, ^ particular interest is the close 'hgreement in pressure near the lead! 

hinge lineTeven thou gli no special shape functions were used** — 





- .4c >11- 2 50 panels k ■ 0..5.1.0 


Figure 26.— Wing Control Surface— Validation Problem Two 
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The^purpose of^P3 is to d^ons^ate ,the use'of RRX fdr*as^^metrical platforms. 
Symmetry conditions employed^n tlxe present study permit the left planform boundary^ 
to be^symmetric,'*anti 8 jrmmetri 5 j.or asymnjjieic. The resolts of an as^metric analysis 
are shown in figures 29 and 30 for a series of skewed •parallelogram planjforms with* 
sweep^^engles ranging from 0° to 60®. The Mach numbee^ds 0.80 and the normalwash 
mode is*that of uniform angle of attack^ Figure 29 ^ows the overall center of pressure 
^fp^susiweep angle, and figure 30 shows^P^e^twn^^lApnt^piof pressure * 





Figure 29.— Overall Center of Pressure for Skewed Wings— VPS 



0 .2 .4 .6 .8 1.0 

ri 


Figure 30. —Sectional Center of Pressure for Skewed Wings— VPS 
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6.6 VALIDATION PROBLEM FOUR' 


L t^r^e, suSso'niffje* twqep^^includi3ig*oi4y the.wpg^exiSl^e'gi4k9b ^ 

I l^sJb^n^calcSlsAed* using ^e^ent metbo^ ^(^tEat of d 

f M^^a^d ^s'^omplir^the result^ o( the preSeA method tl^ gs^of IS>\we f«r 

P '*^the%i»gHitnd%(OfJth| res^Jorifee^Braplitude^with fixed controls for root shear flfrc*, ^ ^ 
^ **^d|^i<ootit>%pd4ng moment, a'iiAviai^p^ccelera’ticA, ^^ectivelv. ^ 
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Figure 32.— Wing Root Bending Moment Frequency Response— VP4 
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Figure 33,-Wingtip Acceleration Frequency Response— VP4 
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6.1 CURRENT STATUS 


f i. 


^ 9 


VI 
» I 


^ 'J -r 

m. f|p^bi|^ity of u^in^ local basis functions for calculating unsteady airloads has»b^n, 

^ ^^e^fully demonstrated. TKe numerical results show 'a degree of smpolhniess ajid, 

^ gccuracy approachijng* that of the global basis function methods while refining ^e, 

^ jjereaUlity an j sm^icitj^ofpahel methods. , ^ 

^ ^or a* gWAi planar ’aerodynanHc''-’surface, Mach .number, reduced frQq>jency,_ ^d_^ 
normalwash distrittition, the computer code developed during the presept ,study will 
predict pressure,* sectional ‘lift, sectional center of ^pressure, overall lift, a/id overall' 
center* of “pi^s^re. If a control is present, values arq calcula^ec^ for tbej>r^^ry surface’ 

* Snd for *the secondary (control)' surface. In addition, force an]^ ^ome^^t ^coeff?ci^t* 

® Sia^rides'can be ealeuleted if desired. Symmetric^, ^^sympietoical^and asymmetrical^ ^ f 

* ^rcfclems can be rum • > • • R 

* ^ ' * * « • » 1 > ♦ » » a 

R Potenttal impsovemente would b%fopnd by qxtening t^e method lo mvAtiple no*np1ai!hiJ • f 


t I 
» I 
• 

-i 

« I 

I 

m I 


6.2 CONCLUSIONS AND RECOMMENDATIONS 


i^f^^aresenj ^copi ^iCi^ (BcSl^lfe^s ’fuhctfciil ^ 

• ljpr{ns^of^c|f(for^i<m Dotenuaf doubled Ian* pi^vfd^ ad d^cA-ate and lasejul^meaiis foiy 4 


» 4i 

nlated iina R | 


At the outset it was recognized that local basis functions similar to the ones noted above 
but formulated in terms of velocity potential doublets, rather than in terms of 
acceleration potential doublets, might provide a more direct growth capability to 
general nonplanar, three-dimensional applications. A parallel steady flow research 
effort by Johnson (ref. 18) utilizing the velocity potential formulation has recently 
demonstrated a superiority over the earlier steady acceleration potential (or vortex 
spline) formulation by Mercer (refs. 11 and 12). However, for unsteady flow, the 
technology for integrating the acceleration potential kernel was readily available, 
whereas the difficulties associated with integration of the velocity potential doublet 
were unknown. Hence, a decision was made at the outset to utilize acceleration 
potential doublets, taking a more conservative approach that would suffice to evaluate 
the salient behavior of local basis functions in unsteady aerodynamics. ••’v 


I 

< 

I 

I 

f 

% 
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I m 
f • 


» * 


I 


6inee«thetative« 1%e upr^b^^n;, o^ qitegr^ting the sul^nic ‘v^oRty p^t^tial *dftibtef 
kemeli has .been a<}(lress^.. This unpublished work, carried out by Johnson, su^e^ts 
integiigtiiig, mostly^n^ closqijl ^orm, the unsteady velocity potential source and doublet 
keqaels yith the, local basis functions used in the work ‘of reference J8. Thus,' it now 
ap^%fs^PQS^ble to work toward the goal hf 'unsteady analysis of < arbitrary 
th]^e^iqiegsu)n^ ^ogfl^rations by implemenffn|; the above mentioaed unsteady, 
kernels into the velocity potential approa^ of ?efSretice 18. • 


I 

I 

• I 

* < 
I 


lii Vi 
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• • 1 t .» » ^4 « I 
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£:©MBUTM?«00»B tfElS<;i« WoJ 

t*' 


« • 


• K • • • C 

!#••••*• * 

I# »«•»••* 

I •••••■ ..T», 

Z 31^ ^Iqw degcrtbeMhe RRX c^e and Ttsxise. 

I « . 4 • 

pi 


. i 

isi'' 

« i 

A.l ORGANIZATION* 


ir ft • • • 

t . ••«««••• *J 

e purpose oC tl^e BR^ (ocj^ i£ to demoiiBt#at4& the Yedsimlity of using the i 

_ii. 2 — ^ t_._i_x^ un8|«e^y. airloads, and tado so in a*manner that could fa^ili^^ 

pn ip th^ fujture to a FLBKSTAB %y^m. ^’he layout of the co^e ^ 
cG^sytewtions in mind^ Figure “A^ depict ftve djst^c; 


the flbw*wfthfti file’ a^u^I *oJ 
: «nits ioi^purpftsdfe ^ planfii^.^ 


i ®p^ts of tfie code, lal^le^d J-\^ whiph ^fipfroximate 

^Ijbi^aratfe t^e woA jjitg rglfi^vj|ly^in^eg(Bi|^e*t «nits tonpurpftseB (Jf planlii^.^ 

I The bulk of the computing lies in calculating the coefficient matrix, denoted by |^C J in 
! part III of figure A-1. From the initial development and numerical analysis viewpoints, 
j part III is the most critical part of the code. Less critical are the input and output 
‘ interfacers, shown as parts II and V. Solution of the linear algebraic equations is shown 
I in part IV. The solutions may be required in terms of the basis function coefficients or 
I in terms of the inverse of the coefficient matrix, and may be either a determinate 
collocated solution or an overdetermined complex least squares solution. 

fcJftA • •OCni'ftilfeNT^ifTRIX AND IJS SOL^TIQN, • > ft * « 

ft *. t ^ ^ * , i aft * « A? 

L^egot^ng by [CnmJ the 'coefffcidht mhtrix appeari^ m equa]uo|^^, ^l)fri 

i_ ;_J i Ai « 1 i^__. 


yidftpc »n4 niisfthe pftessure 


basis function index, we hay^ 


JLBJJi 


( w(x ,y 


T 

Nyy X 1 


-) 


= - -i fc 1 

8tt L mnJ 

N^x Np Np xl 


K) 


(A-l) 


1 

1 


1 


[. Equation A’-l is solved detejjjpingtoly if, and only if, the number of downwgshifxii^Cs | 
j equals the n^mhgr Np of pressure basis functions. The solution may be ove^jj|etavpii*Bd ^ 
in the sdhse of^^mplex least ^uare error by using more downwash*pot^ts th^ b«Bis’* > 
^ function^ ij^bat c^e,4Nyy>Np|and (A-l) i s j~ep l§ced 1^ ^ ^ 


( “ ) - . ji rc^i^rwi [C^i (b„), 

NpXNy NyXi NpXl (A— 2) 
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Figure A- 1.— Program Organization— Preliminary Layout (Schedule Development Phase) 
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The actual calculation of the coefficient matrix is performed as described in section 4.2ji 
However, when the left planform boundary intersects a’ vertical reflection plane (i.e., 
symmetric or antisymmetric, but not asymmetric), the int^ration for Cnin is extended! 
from the support of to its reflection Rn about the left boundary and the value ofj 
Pn or Rp is equal to, or equal to the negative of, the value of P^ on Rn, depending upon 
whether the reflection is symmetric or antisymmetric^ Thuar, for nonasymmetriq^ 
conditions, ^ wa ^ ^ ^ « 


Cmn =ff t - t> dA , 

nm JJ n n pm ’ 


(A- 3) 


n n 

^ jr:^arac{&ristic'<uncti<)n. 




Xn<t) = 


+ 1 
+ 1 
-1 


t e R 


t £ Rj 
t £ R; 


symmetric 

antisymmetric 


(A-4) 


Equations (A-3) and (^^)T^uce t^e numfier of unknowns 'by'a fector»of two when they^ 
apply. The dependence of C mn on both do^wash locationKxm. ynj) an<Lthe Locaii<)n 
the supp^^n o^h^ l^ste function Pn ts deplictfedlnffi^re A^2m * 9 t ft 


r 

I 

I 

I 


✓ 


I 



(Reflective Left Boundary Depicted) 


Figure A-2.—Downwash Point and Pressure Support for Typical Planform 
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_ Q« nyis;^ l^jC jpejfopagd w ^aneC since tire 

I ^cdbfHcMniflnlttaxasi^e#n^4)^l)f Qre^iyej3^i^f^c^o|^ ^i^er, it ?s HeftsSsir^ ft 
I lidAeiniteeKwiiiilti ijj^t^^'^he j&iftl ^t^fti 'fcft 

^d|tiiKi»e<t »eM^Meyay jflMetipqj. . j-gt * 

The; pan^ls*are arra^n^<f iAo^c rdws and 4^8 cnlumnB and are diptipct ftj>n^ ®P®* ^ 
another. The suppo^* h^w^ver, overlap 9n «n essential fashion but«ai^ ^r;^nge4 a^O|| 
into Nc rows and'N^ c^vfinife. ^lie rftimbeS’ of supports •qii^lgtlip Qu^bgr gf ^aij^l^ and^ 
in both cases we numffer theirWrom^l to Np= NcNa starti»g«in4;he first row, proceeding 
left to right, thence to the second row and so on. This numbering system causes the 
index to be given in terms of the row and column numbers by equation (5) for both the 
panels and the basis functions. 


Thus, associated withf 
cqnstituteHhe sg pyiBt 


itk eaclv basis funq^iop ipdex 41 a^sq^ indices*o^aSaefs ftalC ^ | 

^ t]|ay>qgiA f'WciioJt If qp fiflltf 1 I _ * f •. f li 


m * 


N(n) = {(p-1) N + v: (y,v) e A } 
s ri 


(A- 5) 



(Vn) 


(N N(n) = {n: n e TT(n)} 


(A-6) 






f®qua*oHB W)*(A-5fc qp4a(4r6j|hav^bqen, p50gj:aipni^d*iiw RI?k^?rJ^iffe%i?d88i#ed* 
ecdfeoAj/fcrftte «oiftpilitati(# fjf tji^kpmel ^u^ctioi^ ^u^y,^n*tlre ^tBaScfte^tHkyf 
ha^ ikettingspevilige^ ^^u^ ^^e^^l^^t^^l^ ^over* ^ fhe basre TOi^tH^n^ 

auppoidaarq|hfr |)iq§ la 3,acheme Jised in mejrsse i* tgi^ * * * • • • > 




A.4 PANEUNG CONSIDERATIONS 


* • 


" A generic panelisation of a^j ^er ^dynamic sur^c? witli a*co1ntftl Ss ^epictecf in figure • 
A-3. Greater panel density ii®®Oed in regions ^frapid*chSnges*in"pressure. Thus, <* 

along the leadiqg-^g;^, Ji^ |nd trailing-eSge regioiA fclSti^elJ §aAow panels are * 
suggested. „In t^e^case^s^own, wfiere *the ^efft pTamann*b^fcafry«i^ a reflactiire«> 
bounOaiy, ^^atqr ^ensi|^ ^s no^ needed*therO. Ao#e4feif ip£A:antK)llis ipresent^ttlien 
increased {^nel densities are re<^ired^on^bdth%i(!bs%f%hO cAiWoMetfd«gvedeenm^ side 
edges. With onl]^ a few additional intermediate i^iifel honndartospitlK number fff paqpls 
can become significant. ^or*trf^ made exAndle'showw, 99 panels are «s«d- » 

-T 


< 

1 
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Once the planform boundary and panelization have been specified, the type of splined 
pressure functions to be employed are determined in accordance with the edge system 
I indicated in figure A-4. 



Note. 9 chordwise by 
11 spanwise - 99 panels 


Figure A-3.— Possible Panelization for Primary Surface With Control 


f ■ ■■ 

v« 

•Leading Edge 



• Interior Line 

1 


•Control Hinge Line 


• Free Edge 

• Reflective (S or A/S) 

• Interior Line 

• Control Side Edge 



• Control Hinge Line 


• Free Edge 

• Interior Line 
•Control Side Edge 


Figure A-4.— Generalized Pane! Edge System Used in RRX 
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A^PL4N£OHIj(I ^ECmC^fftN 


■•-I 




The planform iS^-spe^ifiMl. aad^a^Knioi;! of a number of suhsurfaceiikaswlepietet^n figure 
A-5. Each ^uBRirface is* a ti:ape^idal ar^ defined by a“nuAber of chord%^s (^gpicted ^ j 
figure ,AJ%. Efdch^on^is defii)ipd by th^ coordinates IJP it^^os^ upstrea®pi pi^nt^nd ^ ' 
length. 'Bpel^caE^nMof th<g ii^pringdiate ch'ord and ^'an Stations as shawn.^ei^uffices ^ 
to d etgini!^ th» pangl^p^gri^ -.1 




Separated surfaces 
not implemented In 
prototype code 


Figure A-5.— Aerodynamical Surface as a Union of Subsurfaces 
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Main Chords (3 shown) 


Column 1 


Column 2 


Intermediate 
chord stations 
(3 shown) 



Figure A-6.— Possible Subsurface Scheme for Prototype Code 
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^ If a discontinuity exists in the spanwise slope of t^e p][an^rm ^u^dary, as depicted in ^ 
! figure A-7, an inadmissible singularity in downj^aslx. results., It is inadmissible in the i 
following sense: the limit of the Kiissner operator^or smooUi edges as the radius of ■ 

' curvature approaches zero doe^ not equal the Kiis^per operator of the limit function (it i 
Ms a discontinuous functionaiiL. Sipce the physical Qow is smoothed across the'boundary ^ 

> layer by the phenomenon of viscosity, it is clear that planfom^reaks’*coul[?f be faired. ^-4 
One such fairing scheme is depicted in figure A-7 in which the fairing width is some 
^fraction, 6p , jjf the wi^th^j^ t^ b<}||.ndeiry tj^p c TOid s. ^ 



Figure A-7.— Fairing Scheme at Planform Breaks 
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: A.6.1 PRESSURE COEFFICIENT 


NP 


ACp(5,n) = Pj^ (5»n) 


(A-7) 
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A.6.2 SECTION LIFT COEFFICIENT 



/ ACp(C,n) 

Cp(n) 

^A^n) 

= ep(n) - |;(nJ / ACpCe.n) de 

5p(n) 


C r 

/ 'SECTIONAL CENTER OF PRESSURE 

I Alfc- » — ' — 




5A<n) 


= c 


5 - tp(n) 


/ ^A^n) -^; r n - y 

” 5p(n) 


^cP 


7^ / 

H r 


CaCh) 


K - ep(n) 

-= ^ AC (e,n) dc 

CA(n)-U(n) P 


Cp(n) 


I A.6.4 SURFACE LIFT COEFFICIENT 

► ___ . » ^ 


= 


■’ll 


ACp(C,n) dA ; A = | | dA 


// 
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(A- 8) 


(A-9) 


(A-10) 


(A- 11) 


(A-12) 
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ACp(C,n) dA ; A = 


II 


dA 


I A.6.5 STREAMWISE CENTER OF PRESSURE ' ^ 


^cP ' c7^ 


1 

“7 

"H 




CACp(C,n) dA 


’cP c 


II 

p 


5ACp(^,n) dA 






A.6.6 SPANWISE'CENTER^F PR15SSURE » * * * * t d 

^ *- ^ I rt i rT- ^ ^*>- * *■ A..,£ f .1 J 


'cP 


' II 


nACp(5»n) dA 


'cP 


II 


nACp(e,n) dA 


"H 


(A-13) 


(A-14) 


(A-15) 


(A-16) 


(A-17) 
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{b } = -8 tt CC ]•' } 

n mn v 

ACp(5,n) = -8itLP^(C,n)J 


(A-18) 

(A-19) 




ACp(5,n) dA 



II <n-nj) 

fj 

Pi 


ACp(C.n) dA 


ACp(5,n) dA 


^ote: f| , m| a nd h&v^^af^c^ eaBfx^,- ^'^,Jac toT^a^ 


(A- 20) 


(A-21) 


(A-22) 
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r 2 . PrbgrairmSd^^uatious^— Siej(^efiTcie«t iStia«es?f "f 1 , ,«an^ 

! to downv^sh: .. ^ - 4 ^ ^ ^ ^ 

i " . - . ^ * 

~ Uic'^l = -«• [|| 


.•antTM/ •witff rebpect 

♦ * «f* 


NPxNW 


‘ I 

NLxNP 


NPxNW 


(A-23) 


[mc«^] = -8. [J| (n-n^ p„<?.n) <ia] [c^]' 


NPxNW 


NLxNP 


NPxNW 


(A-24) 


[aic;„] = -8. [|J p„({,n) dA] [c^]“ 


NPxNW 


NLxNP 


NPxNW 


(A-25) 


' Note; If tt]g least^squac^ ^li^o% isitused (Ny^ > Np), the gjten'efalized invei%e from 
i^us|^bq^u^sti^tq^ fi^ tJip iqdiq^Ltqd ioyecse.ik ^ ^ 

A.6.8 GENERALIZEp ^E^Oipri^lttIQ^I^QRGES« % « m % % ^ ^ 

* ^ ^ ^ ^ ***• ■* ’ 

The generalized aerodynamical^force coefBciqnt delink earlier diaai not been 

programmed Into the RRX For future r^er^nc^ byywQj^ei^the elation to program 
I this importan't' q^ai^it}^ is displayed ^el(]^.^t ^hq)^ld%b^naited^th9it -<the' essential 
practical '’difTicufty 'Yies**in the proper ^nterpgla^on^o^h^stjpcl^ra^ n^>de^ since, in 
I practice, *thfe cOnfeiu'&m*dis^n^^^^^a^^|^d(^ 

*rs ' -*• jf lP„(£,n)J dA CC^r’(^) 

R ® (A-26) 
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